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ABSTRACT 

Precise determinations of the image positions in quad gravitational lenses 
using VLBI can be used to measure the transverse velocity of the lens galaxy 
and the observer. The typical proper motions are //as yr~^, so the time scale to 
measure the motion is ten years. By measuring the dipole of the proper motions 
in an ensemble of lenses we can set limits on the deviation of the inertial frame 
defined by the lenses from that defined by the CMB dipole and estimate the 
Hubble constant. The residual proper motions after subtracting the dipole probe 
the evolution of peculiar velocities with redshift and can be used to estimate the 
density parameter Oq- For N lenses, VLBI measurement accuracies of erg, and a 
baseline of T years, we estimate that the 2-cr limit on the rms peculiar velocity 
of the lens galaxies is 3100((Te/10/ias)(yrs/r)/A^^/^ km s~^, and that the time 
required for the 2-cr limit to reach the level of the local rms peculiar velocity 
Vo,rms is approximately 10A^~^/^(uo,rms/600 km s~^)((Te/10/ias) years. For a ten 
year baseline and = 10 lenses we expect the l-cr limit on the misalignment 
with the CMB dipole to be = 20° or equivalently to obtain an upper limit of 
AHo/Ho < 0.34. 

Subject headings: gravitational lensing - cosmology - peculiar velocities - VLBI 
- astrometry 
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1. Introduction 

Extragalactic proper motions would be a major new cosmological tool if they could be 
measured. For example, an ensemble of distant extragalactic sources defines a rest frame, 
so the average dipole of the proper motions measures our peculiar velocity independent of 
the microwave background dipole (Kogut et al. f993) and local kinematic determinations 
from galaxies or clusters (e.g. Riess, Press & Kirshner f 995; Postman & Lauer f 995; and see 
Dekel f995; Strauss & Willick f995 for reviews). In standard models the directions of the 
two dipoles should agree exactly, and the magnitudes of the two dipoles should agree up to 
the uncertainties in the distances to the sources. If the proper motion dipole and the CMB 
dipole agree in direction, then the alignment falsifies super-horizon sources for the CMB 
dipole (Langlois & Piran 1995; Paczyhski & Piran 1990). Proper motions also limit the 
rotation of the reference frame, where the rotation may refiect a globally rotating universe 
(e.g. Birch 1982). If we assume that the proper motion dipole, and the CMB dipole, 
U0, agree in magnitude and direction, then we can determine Hq by comparing the lens 
redshifts to their distance Hq^ oc D oc Vq/ fi. The residual proper motions, after subtracting 
the dipole due to our motion, measure the peculiar velocities of galaxies at cosmological 
distances. By measuring the proper motions of higher redshift sources, we can begin to 
observe the evolution of proper motions with redshift, thereby having an independent test 
of the predictions of structure formation models and the density parameter Oq- 

Outside the Galaxy, only VLBI has sufficient angular resolution to detect proper 
motions because of the large distances. At a characteristic cosmological (proper motion) 
distance oi D = IDcjpch^l Gpc (1 h'^l Gpc corresponds to z = 0.3 for Oq = 1 and 
Hq = Tbhrs km s~^ Mpc~^) an object with a transverse velocity of io = lO^cuioo km s~^ 
shows a proper motion of // = 0.021ijJiooDqp^ //as yr~^. Since realizable VLBI accuracies 
are ~ 10 //as at cm wavelengths {e.g. Marcaide, Elosegui & Shapiro 1994; Campbell et al. 
1994; Campbell 1995), only velocities ^ 5000-D Gpc km s~^ are detectable in one decade. 
As a result, extragalactic proper motions can be measured only in nearby maser sources 
[e.g. the disk in NGC4258 (Miyoshi et al. 1995) where D = 6.4 Mpc and u ~ 600 km s'^] 
or in superluminal jet sources. Superluminal jets combine both high physical velocities 
[v ~ c) with relativistic aberration (see Zensus & Pearson 1987) to produce proper motions 
between the AGN core and the emitted "blobs" of order mas yr~^. 

Unfortunately, the peculiar and virial velocities of galaxies and clusters of 
100 km s~^ ^ V <^ 1000 km s~^ are too small for direct detection given current life 
expectancies. If the standard error in the position of a source is ag = lOcrgio //as, then the 
time scale for a l-cr measurement is 670cu^QQ-DGpcCeio years. Barring a dramatic increase in 
life expectancy, we need sources with two features missing in typical VLBI sources. First, 
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we need an effect filce refativistic aberration to magnify tlie intrinsic veiocities by a factor 
of 10 to 100. Second, tlie sources must liave nearby reference sources. Proper motions must 
be measured relative to another source, and the more distant the reference source becomes, 
the more difficult it is to control the systematic uncertainties in the relative positions. 

There are VLBI sources that have these features - gravitational lenses. Lenses magnify 
the source, so the apparent velocity of the lensed images is boosted by the magnification 
tensor. Particularly in the quad lenses, the typical magnification of the images is 
several tens, exactly the magnitude needed to make the proper motions measurable on 
"experimental" time scales. The motions in lenses include not only that of the observer 
and the source, but also the lens galaxy. In fact, the motion of the lens galaxy probably 
dominates over the source motion. The lens galaxy lies at fairly low redshifts 0.1 <^ z <^ 0.5 
(some lens galaxies can be very close to the observer - the lowest is 2; ~ 0.031 for the 
optical quad lens 2237+0305, Huchra et al. 1985), so they are both closer to the observer 
and have larger peculiar velocities than the sources. Finally, by measuring proper motions 
in multiply imaged sources, we can measure the relative positions of the images rather 
than their absolute positions. This eliminates the need for an external reference source, 
and minimizes the systematic errors because all the sources/images are confined to a small 
angular region. Moreover, by measuring the relative motions of the images, we typically 
double the apparent velocity because pairs of highly magnified images generally have 
anti-parallel motions. If we scale from existing radio lens surveys (MG, Lawrence et al. 
1984; Burke 1990; JVAS, Patnaik et al. 1992a; and CLASS, Myers et al. 1995) we expect 
20 fiat spectrum lenses brighter than 100 mJy, 40 brighter than 50 mJy, and 300 brighter 
than 10 mJy on the sky, and the bright samples are evenly divided between two and four 
image systems. Proper motions in gravitational lenses have previously been discussed 
in the context of superluminal sources. Chitre & Narlikar (1979, 1980, also Chitre & 
Saslaw 1989) discuss gravitational lensing as a possible cause of superluminal motions, and 
Gopal-Krishna & Subramanian (1991, 1996) discuss it as a means of proving that compact 
double radio sources are gravitational lenses. 

In §2 we discuss the CMB dipole, our peculiar velocity, and the dependence of peculiar 
velocities on redshift and cosmology. In §3 we discuss the magnification of proper motions 
by gravitational lenses. In §4 we discuss the current VLBI accuracies and the feasibility of 
measuring the magnified proper motions. In §5 we discuss using the mean proper motions 
of lensed VLBI sources as an inertial reference frame, and in §6 we discuss using the 
residual proper motions to study the evolution of peculiar velocities and discuss systematic 
biases in the measurements. In §7 we consider effects other than peculiar velocities that 
may contribute to the observed magnified proper motion. Finally, in §8 we summarize our 
results. 
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2. Peculiar Velocities 

The advantage of peculiar velocities as tools for studying large scale structure is that 
galaxies are used only as "test particles" moving in the gravitational field, thereby avoiding 
the unsolved problems associated with galaxy formation and the type of matter inducing 
the acceleration. Existing techniques determine the solar system peculiar velocity from the 
CMB dipole anisotropy to be Vq = 386.6 km s~^ (Kogut et al. 1993), and the peculiar 
velocity of the Local Group {e.g. Yahil, Tamann & Sandage 1977) to be vlg — 600 km s~^. 

The determination of peculiar velocities is strongly limited by the accuracy of distance 
indicators. A typical distance indicator has a fractional error e, where the value of e ranges 
from 0.05 for Type la supernovae (Riess et al. 1995) to 0.21 for the D^-cr relation of 
elliptical galaxies (Dressier et al. 1987). On top of this intrinsic scatter in the correlation 
between the distance dependent and the distance independent quantities, there are also 
systematic errors that can only be removed statistically. For the measurement accuracy to 
be comparable to the radial peculiar velocity, v^msl V^-, we must observe 

, ecz \ f e f cz \^ ( 600 km s~^ 

iV ~ 3 ~ 19 



0.15/ VlOOOOkms-i 

independent objects at redshift z. This exceeds the local number of IRAS galaxies in a 
(5 /ifs^Mpc)"^ volume by a factor of ~ 2.5. Even at distances much smaller than 10000 km s~^, 
the intrinsic scatter is large, and averages over large volumes (~ (10 /ifs^Mpc)"^) are used to 
reduce the noise level. As a result we measure only coarse averages of the peculiar velocity 
over many galaxies rather than the peculiar velocities of individual galaxies. We estimate 
the raw, unaveraged peculiar velocity from the residuals about the mean peculiar velocities 
corrected for the estimated intrinsic scatter in the Hubble diagram, or from the pair-wise 
velocities of galaxies. The best observational estimates for the rms peculiar velocity range 
from 430 km s~^ to 700 km s~^, although estimates can be as low as 300 km s~^ and as 
high as 1800 km s~^ depending on the sample and the analysis method (see Strauss, Cen & 
Ostriker 1993). 

Structure formation theories predict the evolution of the rms peculiar velocity, and 
when coupled with the COBE normalization (Gorski 1994, Gorski et al. 1994, White & 
Bunn 1995) for the power spectrum, predict the amplitude of the rms peculiar velocity 
today. Linear theory (e.g., Peebles 1993, §21) predicts that 

^oV^(»o,Ao) |- 
2 7r^ Jo 

where a zero subscript denotes a value determined at the current epoch, Hq is the Hubble 
constant in units of km s~^ Mpc~^, P{k) is the COBE normalized power spectrum, and 
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the function / (Lahav et al. 1991) is 
/(Oo,Ao)^Or + ^(l + y) and /(Oo,Ao,^) 



Oo(l + ^) 



Oo(l + ^)3-Oofc(l + ^)2 + Ao 



0.6 



(3) 

where the second equation neglects the second term of the first, and where 0^ = 1 — — A 
is the "curvature density" of the model. For example, in standard COBE normalized CDM, 
the predicted rms peculiar velocity today is Vrms,o = 760 km s~^ with a 100 kpc filter. The 
rms peculiar velocity evolves with redshift as 

_ /(Oo,Ao,z) 1 

/(iio,Aoj Vl + -2 

and the rms transverse velocity causing the proper motions is just cOrms = (2/3)^/^Urms(-2)- 
We must choose whether to use the theoretical models to predict both the evolution and 
amplitude of the velocities (eqns. (2), (3) and (4)), or simply to evolve the locally observed 
peculiar velocities to high redshift (eqn. (3) and (4)). 

The rms proper motion as seen by an observer in the CMB rest frame is //cmb = ^/-Do5, 
where Dqs is the proper motion distance between the observer and the source. We use 
distances defined by 



{lOofcr/' [(1 + zfil + noz) - z{2 + ^)Ao]-'/' dz} 



r smn 

i)(0o, Ao, ^1, z,) = — ^' ' 1^ I,,; ' ^ ' ' (5) 

-no |"ofc| / 

where the function sinn(a;) is sinh(a;) for an open universe (0^ > 0), sin(a;) for a closed 
universe (0^ < 0) and simply x for a fiat (0^ = 0) universe (see Carroll, Press & Turner 
1992). When z\ = ^ this is the proper motion distance to redshift z^} 

Given the apparent angular velocity in the CMB frame, the effective proper motion of 
a galaxy at z due to the combined solar system CMB dipole and the source velocity cu^ is 
given by 

Vq-{z-Vq)Z CO, 
/^e = J, J^. (6) 

J^os J^os 

If, however, we are interested in the relative proper motion of a lens galaxy, L, at redshift 2;;, 
and a source galaxy, S, at redshift Zs > zi, as observed by an observer, 0, then the effective 



"'^For zi 7^ it is not the proper motion distance of an object seen at Z2 by an observer at zi - it lacks a 
redshift factor. The distances have the convenient algebraic property for aflat universe that D13 = -Di2+-D23- 
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relative proper motion becomes 

= 7^ ( ^' - 7^^« J (V 

J^OL J^OS J^OL V L>os / 

(Kayser, Refsdal & Stabell 1986, Miralda-Escude 1991) where uji and ujs are the transverse 
velocities of the lens and the source. For zi < 0.3 and Zs — 2 we obtain in the range 
0.1 < Oo < 1 

— 77^ < 0.20 (8) 
^/ L>os 

so we can usually neglect the source galaxy motion in eq. (7) and simply use the proper 
motion given by 

l^'e = -^Y^@-{z-Vq)z\-^ . (9) 

UoL y Uos > 

This approximation, however, is not valid in the presence of large internal motions within 
the source (see §7). As discussed in the introduction, these proper motions are of the order 
of ~ 0.1 //as yr~^ and thus can be observed only if appreciably magnified. 



3. Lens Magnification of Proper Motions 

In this section we first explore simple, generic models for the expected magnification of 
proper motions for fold caustics, cusp caustics, and symmetric lenses, and then we examine 
simple models for several of the known VLBI lenses. The general lens equation is 

tT=f-VV', (10) 

where x and u are the angular coordinates of the lens plane and the source plane 
respectively, and %l) is the two-dimensional effective potential of the lensing mass distribution 
(see Schneider, Ehlers & Falco 1992). The inverse magnification tensor is the Jacobian of 
the lens mapping 

M^^ = %-V.,; V-.-^^, (11) 

and the magnification diverges on the critical curves where the determinant of the Jacobian 
is zero, Af~^ = lAf^^^l = 0. Taking /Te of eq. (7) as the effective proper motion of the lens, 
the proper motion /T/ of an image located at xi becomes 

/Tl=M(fj)/re (12) 

where M is the magnification tensor with components Mi,j. 



-7- 



When we need a concrete example, we use a singular isothermal sphere in an external shear 
field, for which 

tjj = br -^-fr^ cos 2(9 - 9^), (13) 

where b = 4:Tr[cr / cy D^s / Dqs is determined by the velocity dispersion of the dark matter 
cr (cr ~ 220 km s~^ for an galaxy), and the shear is characterized by a dimensionless 
strength 7 (typically 0.05 < 7 < 0.15) and an angle 9^. This potential is a reasonable, but 
not perfect, model for lens statistics and lens models (see Kochanek 1991, 1996). 



3.1. Motions In Particular Lenses 



A source on one side of a fold caustic has two images on either side of the critical curve 
separated by 



\Ax\ = 2 



, =-i\ ni/2 
2tr M Am 



|e • VM-i| 



(14) 



where VAf ^ is the gradient of the inverse magnification on the critical curve, e is the 
singular eigendirection of M on the critical curve, and Am is the distance of the source 
from the fold caustic. Two images are present for Am > 0, and none are present for Am < 0. 
To first order in Am, only the component of the source velocity perpendicular to the caustic 
[du±/dt) produces changes in the image separations. For the simple model potential (13), 
the time derivative of the distance between the images is 



dAx 



dt 



46 



1 



Ax 37 sin 2^ 



du± 



dt 



1 +7 cos 26* + 0(7^ 



(15) 



to first order in the shear 7 for the model potential (13). Ignoring the small multiplicative 
corrections for the shear, the perpendicular velocity is magnified by (46)/(3Aa;7 sin 2^). 
The fold approximation is invalid for sin 2^ = 0, where the caustic has cusps. 

The other generic feature of lens caustics is the cusp, where three rather than two 
images merge. We assume that the cusp is axisymmetric with respect to the Ui axis, that 
its tip is at the coordinate origin, and that it points towards Mi — )• 00. A source on the axis 
with Ml < has three images. One image is formed on the Xi axis and has negative parity, 
the two other images have positive parity and are located symmetrically on opposite sides 
of the Xi axis. The relative motion of the two images with positive parity for a source on 
the symmetry axis is (e.g. Schneider, Ehlers, & Falco 1992) 



(16) 



d\Ax\ 




#122 


1 


dui 


dt 




V'122 + |V'2222(1 - V'll) 


|Aa;| 


dt 
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This can be evaluated for any lens model, but we can find a more general solution for 
our simple potential (eq. (f3)) that includes this limit by explicitly solving 

for the image positions when the source is on the Ui axis for 6^ = 0. If < b 
then the lens produces at least two images located at Xi = {{ui + b)/[l — 7),0} and 
X2 = {{ui — b)/[l — 7),0}. If < Uc = 276/(1 + 7) the lens equation has two additional 
solutions at X3 = {— m/27, j/d and X4 = {— m/27, — j/c} where j/c = b[l — m^/m^)"^/^/(I + 7). 
This case is illustrated in Figure I. The relative velocity of images 3 and 4 is 

dAx34 _ 1 2 b 
dt (I + 7) 7 Aa;34 

where the cusp limit is AX34 <^ b. The factor under the square root suppresses the relative 
motion of images 3 and 4 in the limit of a symmetric quad lens {A.X34 ~ 26/(1 + 7), u <^ Uc), 
where the source is near the origin and the images are arranged in a cross with images 3 
and 4 forming the top and bottom images (as in 2237+0305.) If we look at the relative 
velocity of images 2 and 3 for u <^ then 

dAx23 1 
dt ~ 27 

so that for the symmetric lenses there is a large magnification of the relative motions of the 
images on different axes. Comparing these equations, and ignoring corrections of order 7, 
the typical magnification of proper motions in a four image lens is of order b/2^fAx ^ 7~^. 
Figure 1 shows a time sequence for a typical lens at zi = 0.3 moving at 500 km s~^ along 
the symmetry axis of the potential t/j defined by 6 = 1.0 arcsec and 7 = 0.1. Notice the 
very rapid motions of images 3 and 4 compared to images 1 and 2. In the four-image lenses 
we would measure the highly magnified relative motions of 3 and 4, while in the two image 
systems we can only measure the much slower motions of images 1 and 2. ^ 

To illustrate these analytic estimates, we took the observational data for four VLBI 
lenses (MG 0414+0534, Hewitt et al. 1992; B 1422+231, Patnaik et al. 1992b; and CLASS 
1608+656, Myers et al. 1995), and fit them with our simple potential model. We do not get 
perfect fits to both the relative positions and fiuxes of the images, but we do get reasonable 
quantitative estimates of the sensitivity of the lenses to proper motions. From the fit we 
determine the magnification tensor at the position of image /, Af(x/). If the effective source 
position relative to the lens shifts by u, then the image position shifts by Sxj = M{xi)8u. 
We can measure either the relative motions of the images Sxj — Sxj or the absolute motions 
of the images Sxj. Figure 2 shows the polar diagrams of motion sensitivity for the three 



'l+-f)Ax34 

2b 



dui 



dt 



:i7) 



dui 



dt 



18) 



^In fact, in this model, images 1 and 2 show no relative proper motions! 




source plane image plane 

Fig. 1. — Proper motions. The left diagram shows the source plane, and the right shows 
the lens plane. The scale of the source plane is magnified by a factor of two. The dashed 
lines in the source plane show the outer radial pseudo-caustic and the inner astroid caustic, 
and the dashed line in the image plane shows the tangential critical line. The source and 
image motions are shown by the solid lines, and the points are spaced at time intervals of 
10^ years for a lens at z = 0.3 moving at 500 km s~^ and a lens critical radius of 6 = 1.0 
arcsec. It takes 9.6 X 10^ years to cross the critical radius. The images start at the image 
number labels for a source at the origin. 

known four-image VLBI lenses, and the classic double lens 0957-|-561. We show both the 
maximum absolute motion of any single image, and the maximum relative motion of any 
image pairs as a function of the direction of the effective velocity. Note the different angular 
scales in the different panels. 



3.2. Statistics of Proper Motion Magnification 

For the statistical analyses of §5 and §6 we need expressions for the mean square 
magnification of proper motions. In the isothermal model (eqn. (13)) the eigenvalues of 
the magnification tensor are effectively Ai = 1 and A2 = Af/ where Af/ is the magnification 
of image /, so if we know the magnification probability distribution, we can estimate the 
average magnification of proper motions. If we assume that the total magnification Af ^ 1, 
then we can allow Ai = 0. For the four image lenses, the mean total magnification of all four 
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Fig. 2.— Polar sensitivity plots for proper motions in MG 0414+0534, B 1422+231, CLASS 
1608+656, and 0957+561. The angle is the direction of the effective source velocity, and 
the magnitude is the resulting proper motion scaled to an effective velocity of 600/^75 km s~^ 
and the listed lens redshift. The solid curve shows the maximum relative proper motion 
between the images, and the dashed curve shows the maximum absolute proper motion of 
any image. The solid vectors show the direction and magnitude of the effective velocity 
and proper motions corresponding to the projection of the CMB dipole velocity on the lens 
(north to top, east to left). The heavy solid circles show the unmagnified proper motions 
expected for a 600 km s~^ velocity at the lens redshift. Note the different angular scales in 
the different panels. 
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images is (Af) = Cj^ where C = 3.6, and the integrated probabihty distribution for the 
total magnification is well approximated by P(> M) = ((M)/2M)^ with {M)/2 < M < oo. 

For high magnifications, the image proper motions are dominated by folds, so the 
effective magnification of the proper motions is ~ M because two images dominate the 
magnification with Ma — Mj, ~ M/2 and Mc — <^ M and the magnification tensors of 
the two highly magnified images will be nearly aligned. Thus the mean magnification of the 
proper motions is 

/■°° dP C 

(^-) = / ^^7^^ = -' (19) 
J{M)/2 dM 7 

and the mean square magnification is 

1/2 

(20) 

The latter integral is logarithmically divergent in the maximum magnification cutoff of 
Mraax ~ bjOg ~ 10"^(mas/^s) where h is the critical radius and 6^ is the source size. For 
statistical calculations we will use the approximation that the mean square value of the 
magnification eigenvalues that enters the statistical calculations is 

{\l + \l)^{M,Y = ^. (21) 

This is an underestimate of the true proper motion magnification because we neglect 
the logarithmic divergence, multiple image pairs, and magnification bias. Magnification 
bias makes fiux-limited samples have higher mean magnifications than the mean of the 
magnification probability distribution. 



{Ml) 



Mrriax dP 26 

(M>/2 dM ^ 9,s{M) 



4. VLBI Measurements 

The observational issue in this project is the accuracy with which image positions can 
be determined using VLBL We define ag to be the l-cr error in one coordinate component, 
and scale all estimates to ag = 10 //as. The simplest measurement of the proper motions 
consists of two epochs separated by time T with rms position errors erg, and this leads to 
rms measurement errors for a component of the relative proper motion of one image pair of 

Additional measurements at intermediate epochs, and averaging over the independent 
image baselines will reduce the uncertainties, but we scale our estimates using this simple 
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measurement model. The estimates we make in §5 and §6 will overestimate the uncertainties 
for any given value of erg. 

Few astrometric VLBI observations closely resemble the problem of determining 
relative proper motions in gravitational lenses, because in most astrometric studies the 
errors are dominated by effects arising from the separation of the source and the nearest 
reference source. For example, Guirado et al. (1995) estimated the separation of 1928+738 
and 2007+777 (separated by 5 degrees) to an rms accuracy of 200-300 //as, where the 
uncertainties were dominated by the tropospheric model and the position of the reference 
source. Marcaide, Flosegui & Shapiro (1994) estimated the relative proper motion of 
1038+528 A and B (separated by 33 arcseconds). In this case both sources can be observed 
simultaneously, as is the case for the lens sources. Such sources can be "phase connected" 
and the possible precision increases greatly. The statistical uncertainties in the source 
positions were ~ 2 //as at 3.6 cm and ~ 10 //as at 13 cm, leading to a nominal proper motion 
of 31 ± 2 //as yr~^ between 1981 and 1983. Unfortunately, an analysis of systematic errors 
in the measurement suggests that the true errors are nearly ten times larger, and their final 
estimate of the uncertainty in the proper motion was ±22 //as yr~^ at 3.6 cm. Only one 
lens, 0957+561, has accurate, multiple epoch, astrometric data (Campbell 1995, Campbell 
et al. 1994). In this system (see §3) there are two images separated by 6 arcseconds, so it is 
even easier than in the case of 1038+528 A/B to determine the relative positions. At 6 cm 
the typical statistical uncertainty in the relative positions was 13 //as for each coordinate, 
but the actual scatter in the data (including earlier results by Gorenstein et al. 1984, 1988) 
was two to three times larger. 

If we reject the possibility that the apparent motions in 1038+528 A/B and 0957+561 
are physical (they correspond to motions with ^ c, but see §7), then current measurement 
accuracies are of order 20-30 //as. Fortunately, the quad lens systems are probably less 
sensitive to systematic errors than either of these systems, because the image separations 
are smaller (2 arcsec or less). This reduces the magnitude of systematic errors that are 
first (second) order in the separations by another factor of 3 (9) relative to 0957+561. The 
presence of four images allows the systematic errors to be estimated from even two epochs 
of data, because the best fit proper motion specifies only two of six relative motions. The 
residuals from the best fit proper motion are a combination of the intrinsic errors in the 
VLBI measurements and errors in the magnification tensor estimates. Space- VLBI and the 
use of higher observing frequencies may also reduce the uncertainties by improving the 
intrinsic resolution of the observations. 
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5. The Mean Proper Motion and The Inertial Reference Frame 

For each lens we measure the changes in the relative image positions over a time 
baseline T. In a lens with n = 4 images, there are n[n — l)/2 = 6 independent relative 

positions. The change in the relative positions of images / and J is T Afijj, where the 

relative proper motion is A fijj = jli — jlj. We assume that the magnification tensor at each 
image M[x) is exactly determined by a lens model, allowing us to solve for the effective 
proper motion jll in each of / = 1 • • • lenses. 

To estimate the expected signal-to-noise ratio for detecting the proper motions we 
consider the case of a single image pair. The magnification of the relative proper motion 
is determined by the eigenvalues Ai and A2 of the difference in their magnification tensors 
(Af/ — Af j), and the relative orientations of the magnification eigenvectors and the peculiar 
velocity. If we average over the orientation of the lens, the location of the lens, the position 
of the source, the lens peculiar velocities, and use the ensemble average {Xf + A2) = C'^^~'^ 
of the proper motion magnification (see §3.2) , then the typical signal-to-noise ratio for a 
single lens is 




0.0031 vq 



{\1 + \IYI' veT 
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(23) 
(24) 



In an Oo = 1 universe with 7 = 0.075, zi = 0.3, = 2.0, Vj-msfi = 600 km s~^, and 
Vq = 369 km s~^ we obtain an rms signal to noise ratio of ~ 0.25/i75(10/ias/(Te)(r/yrs). A 
temporal baseline of T = 12/^^5^ (crg/lO/xas) years leads to a 3-(T detection for a typical quad 
lens. 

In each lens the proper motion is the sum of a proper motion dipole produced by our 
motion relative to the inertial frame and the randomly oriented proper motions of the lens 
and the source. We can determine our proper motion by fitting a model proper motion 
dipole il[^{zi,A) to the measured proper motions in an ensemble of gravitational lenses, 
where A = (Oq, Aq, Hq^ t^o,rms, -P(^)) contains all the cosmological parameters. Formally, 
we would fit the data using a statistic 



X 
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/=i 
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pec m 



{zi ; A) 



(25) 



The additional term of fipecm denominator includes the noise in determining the mean 

dipole from the rms proper motions of the lens, and the e; are the measurement uncertainties 
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translated to the uncertainties in the unmagnified proper motions (e; = e;/(|Af|)). Neglecting 
the noise Irom the lens proper motions, the mean square value of the is simply N times 
the first term of equation (23). 

Two simple results are the uncertainties in estimating the direction of the dipole 
velocity and the value of the Hubble constant assuming that the CMB dipole and lens 
dipole are identical. If we define as the uncertainty in the angle between the true dipole 
velocity and the model dipole velocity, and AHq/Hq as the fractional uncertainty in the 
value of the Hubble constant, then the l-cr limits (Ax^ = 1) foi" ^ lenses with the noise 
dominated by the measurement error are approximately 



Ai^o . , , [Ax'Y' f 7 \ f '^e \ f Dos \ /10yrs\ 



or sm . 



where {D\g j D^qj^D^q^) = G/Dq^ is an average over the lensing cross sections. For = 10 
lenses the l-cr limit is A^ = 20° and the 2-cr limit is A^ = 43° after ten years. This estimate 
is a lower limit for a given value of erg, because the uncertainties for a particular lens can be 
reduced by considering all n[n — l)/2 relative image positions and by making more than two 
measurements of the image positions, and furthermore we are underestimating the typical 
magnifications. Longer observing time scales do not help because the error is eventually 
dominated by the real proper motions of the lenses rather than observational uncertainties. 
The ultimate limit on the uncertainty of the proper motion dipole is constrained by the 
total number of lenses. As can be seen in Figure 2, the proper motion due to the peculiar 
velocity of the lens is typically larger than the proper motion due to the peculiar velocity of 
the observer. 

Figure 3 shows the estimated dipole direction using the location of known VLBI lenses. 
We drew the peculiar velocities from a Gaussian distribution with Vrms,o = 600 km s~^ 
scaled to zi = 0.3 in an = 1 universe, used average magnification of C/7\/2 = 34 with 
C = 3.6 and 7 = 0.075, and Gaussian error distribution with ag = 10 //as yr~^ over a time 
baseline of 10 years. We ran 50 Monte-Carlo realizations and determined the mean dipole 
direction and its dispersion. We also show, on the same plot, the inertial frame dipole as 
obtained if all lenses rotate about the Galactic pole with rotation velocity equal to their 
expected rms two-dimensional peculiar velocity (thinner cross). 



6. Peculiar Velocities at Cosmological Distances 



Once the observer's peculiar motion is subtracted, the remaining motions are due to 
the peculiar velocities of the lens and of the source. For an Oq = 1 universe, the peculiar 
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+90° +80° 




Fig. 3. — Predicted inertial frame dipole direction (bold cross) as calculated by 50 Monte- 
Carlo runs on known VLBI sources locations with varying velocities and measurement errors. 
All lenses are considered to be at 2; = 0.3. The heliocentric CMB dipole is marked by the 
filled dot. A thinner cross is plotted for comparison, when we make all the lenses rotate 
about the Galactic pole with (uf{z)Y/''^. The figure is in Galactic coordinates. The dotted 
line marks the celestial equator. 



(see eq. 
^^^lrms9^{^)lDls where 



velocities scale with redshift as {ijj'^{z)) = ^orms/(l + where cUq^^^ = (2/3)uq, 
(3)), so the lens cross section weighted average of {cof / Dqj^) 
g[z) ~ (1 + ZsY^^ for Oo = 1- The contribution of the peculiar velocity of the source galaxy 
is negligible compared to that of the lens (see eq. (8)), because the relative contributions 
are 



0.1 1 + 



-11/9 



(27) 



with the correct cross section weighted averages of the velocities. Thus the upper limit on 
the mean square peculiar motions scales as 
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and for N lenses the time required for the 2-a upper hmit (Ax^ = 4) to reach the level of 



r = 10 
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(29) 



Figure 4 shows the expected rms proper motions over ten years in a specific, simplified, 
family of large-scale structure models. The lenses are all at a redshift of zi = 0.3, and the 
proper motions were all magnified by the typical value of 34. We used a COBE normalized, 
CDM power spectrum for open cosmologies ("OCDM" models with Aq = 0) and assumed 
that the lenses were a fair sample of galaxies. Models with small Hubble constants and 
low matter densities have very small rms proper motions (an rms of 1.4 //as after ten years 
for Oo = 0.2, Hq = 40 km s~^ Mpc~^) while models with high Hubble constants and high 
matter densities show large proper motions (an rms of 150 //as after ten years for Oq = 1-0, 
Hq = 80 km s~^ Mpc~^). Small Hubble constants reduce the proper motions because the 
distances to the lenses are larger, while low matter densities produce smaller rms peculiar 
velocities. Similar result holds for a fiat universe model (ACDM) with Oq and Hq as 
variables. In these models we assume the theoretically calculated rms peculiar velocities 
from eq. (2), not the observational estimates of the rms peculiar velocity locally. 

We do not escape systematic problems by using lens peculiar velocities. In particular, 
the lenses are a biased sample, because most lenses [roughly 7 out of 8, Kochanek (1996)] 
are formed by E and SO galaxies rather than spirals - the higher masses of the E/SO 
galaxies more than make up for their lower number density compared to spiral galaxies. 
The morphology-density relation (Postman & Geller 1984, Dressier 1980) shows that E/SO 
galaxies are ten times more likely to be in or near clusters of galaxies, so the sample 
defined by the lenses is biased to represent the pattern of velocities in denser regions. 
Since the pair-wise velocity dispersion of galaxies increases with the density of the sample 
environment (Marzke 1995), the rms peculiar velocity of lens galaxies will be higher than 
average. This occurs even though the higher density environment has only weak effects on 
the lensing properties of the galaxy - we know from the morphology of the lenses that the 
lens galaxies do not lie in the inner ~ 1 h^^Mpc of any cluster. However, the velocities 
will be strongly modified if the galaxy lies within 3-5 h^^Mpc of the cluster center."^ Thus 
the peculiar velocities of lens galaxies will typically be higher than the average peculiar 



^The perturbation a cluster causes on the lens model depends on the projected gravitational acceleration 
produced by the cluster at the radius of the lens galaxy and is related to the weak shear of the cluster. 
Outside the cluster core, this is only a small perturbation on the lens model. The perturbation in the 
velocity of the lens galaxy by the cluster is much larger because the velocity perturbation is the time integral 
of the weak acceleration. 
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Fig. 4. — Expected measured piopei motion deviations about the mean (dipole). Velocities 
and proper motion distances are lor zi = 0.3 in an open universe CDM (OCDM) model with 
the average magnification ol ~ 34. Contours are spaced Irom the minimum at the lower left 
corner (1.27/ias) by A/i = lO/xas . 



velocities of randomly selected galaxies or field galaxies. This makes measuring the proper 
motions easier, but their interpretation more difficult. 



7. Other Sources of Proper Motions 



Peculiar velocities are not the only possible source of proper motions in lensed VLBI 
sources. Proper motions may also be produced by motion in the source AGN, and possibly 
by microlensing in the lens galaxy. Although the peculiar velocity of the galaxy containing 
the AGN is negligible, the velocity of the source may be dominated by internal motions 
near the AGN. If typical sources have effective internal three-dimensional velocities u^, 
then the internal motions dominate the proper motion residuals in an Oq = 1 model if 
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Vs ^ (l%)^^^^o,rms = 3.9(1 + ^;s)^^^^o,rms- H wc take iOQ^rms = ^^2/3 600 = 490 km s ^ then 
source motions dominate if they are larger than Vg ^ 1700 km s~^. 

The only effect that mimics peculiar velocities over long time scales is motion of the 
supermassive black hole powering the AGN relative to the source galaxy. Objects in orbits 
smaller than 3 X 10^'' cm at a distance of 4 Gpc have maximum proper motions smaller 
than 10 //as and have no effect. Orbits larger than lO^^Afs cm, where the mass producing 
the orbit is Af = lO^AfsAf©, have source velocities that are only 10% corrections to the 
expected peculiar motions. Between these two limits, only the wider orbits are physically 
plausible either because of limits from gravitational radiation for black hole binaries, or 
because of limits on plausible stellar densities for a black hole orbiting in a central star 
cluster. The maximum plausible orbital effects are comparable to the effects of the lens 
peculiar velocity, but in most cases the contamination of the peculiar velocity from motion 
of the black hole powering the AGN should be negligible. 

Apparent peculiar velocities can also be produced by the motion of the emitting 
material relative to the black hole. The motion may either be orbital motion or non- 
gravitational motions produced in the jet. Radio emission from orbiting material is unlikely 
to be important. If we suppose that the emission is created by N patches with velocities 
generated by the gravity of the black hole, then the maximum characteristic velocity 
dispersion of the mean centroid of the emission is of order 

a, ~ 3000N-'/'Ml''D-ali' km s'^ (30) 

where (Ts is the angular size of the emission region. Unless ~ 1, we do not expect the 
centroid of the radio emission region to show a large enough velocity relative to the AGN 
to be important, and there is no reason for orbiting material to reach this limit. 

The remaining case, and the one known to be important from the existence of 
superluminal sources, is motion of material in the AGN jet. The effects of source motions 
and the possible relations between gravitational lensing and normal superluminal motions 
are discussed by Chitre & Narlikar (1979, 1980, also Chitre & Saslaw 1989) and Gopal- 
Krishna & Subramanian (1991, 1996). For our purposes, we need not worry about classical 
superluminal motion. A lensed classical superluminal source would be easily detected, since 
the magnified motion of the superluminal components could be tens of milliarcseconds per 
year. Suppose the core AGN has fiux fi and a blob is emitted with fiux /2 and apparent 
transverse velocity of io. Then the apparent motion of an unresolved source is iof2/{fi + /2)- 
Sufficiently fast and bright components are detectable as normal superluminal motion, so 
the primary problems are superluminal but very low fiux components, and "sub-luminal" 
bright components. Since such motions are, by definition, not directly observed, the only 
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way to eliminate them as a source of error is to monitor the source for extended periods of 
time - the true proper motions will continue to grow, while internal source motions will 
be an extra source of stochastic noise. If the motions seen in 1038+258 A/B or 0957+561 
discussed in §4 are real rather than systematic errors, then proper motions will be trivially 
detected in lenses. A quad lens should magnify the 10 — 20/ias yr~^ motions to 300 — 600/ias 
yr-^ 

The other process that may affect the proper motion measurement is microlensing of the 
source by the granularity in the potential of the lens. Gopal-Krishna & Subramanian (1991, 
1996) have suggested that microlensing may produce observable effects in compact double 
radio sources. The overall magnification of the source is produced by a combination of the 
smooth potential creating the multiple images, and microlensing by the fiuctuations in the 
gravitational field from the stars in the lens galaxy. The effects of motion have mainly been 
discussed in the context of microlensing the optical cores of AGN (Kundic & Wambsganss 
1993; Wambsganss & Kundic 1995; Kundic, Witt & Chang 1993). While microlensing is 
detected in several optical lenses as time variations of the optical fiux (e.g., Irwin et al. 
1989; Corrigan et al. 1991; Racine 1992), it is impossible to determine the velocity of 
the images from the light curve with any precision even with perfect sampling. The fiux 
variations are a combination of the velocity, the mass spectrum of the lenses, the size of the 
source, the surface matter density and the shear. It is also a stochastic process, so unlike 
the proper motions of the primary images there is no signal that monotonically increases 
with time. Microlensing by solar-mass stars acts on a characteristic angular scale of 10~^ 
arcsec while the size of typical VLBI sources is of order 10""^ arcsec. Hence microlensing by 
individual stars affects a fraction of 10~^ of the source size and produces a negligible effect 
on the total fiux of the sources. Although the combined effect of an ensemble of stars can 
affect a much larger region of the source (typically about ten times larger, Katz, Balbus, 
& Paczyhski 1986), the magnification of the total fiux by microlensing remains negligible. 
In this case, microlensing cannot induce any detectable motions in the images. For the 
characteristic microlensing scale to be of order 10""^ arcsec, the microlenses would have 
to have ~ 10^ Mq, typical for globular clusters, but the optical depth for microlensing by 
globular clusters is very small. Even then, the random motion of the microlenses within 
the lensing galaxy is of order 100 km s~^ or about an order of magnitude smaller than 
the transverse velocity of the macrolens, rendering possible microlensing-induced image 
velocities very small. Although the velocity of microlensing caustics can be much larger 
than that of the microlenses, random velocities of microlenses in the lensing galaxy are 
on average only ~ 25% more efficient in producing time variable magnifications of the 
source than the transverse velocity of the galaxy as a whole (Kundic & Wambsganss 1993). 
Finally, microlensing would only affect the fiux of one source image at a time so that the 
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effect of microfensing coufd be identified by using the reiative veiocities of aii source images. 



8. Conclusion and Discussion 

Detectabie proper motions in tlie quad ienses observabie witli VLBI are inevitabie, 
and tlie oniy issues are tlie time scaie required to detect tlie motions, and wlietlier tlie 
motions wiii be dominated by tlie peculiar velocities of the observer and the lens galaxy or 
motions in the source AGN itself. Although we focused on the study of peculiar velocities, 
either result is physically interesting. Useful limits can be set given temporal baselines of 
10 years, and it is actuarially certain that proper motions will be measured in all of the 
VLBI quads before the authors of the paper retire. For the moment, the only observational 
requirement is that VLBI observations of the quad lenses should be made keeping the goal 
of determining accurate proper motions in mind. The proper motions expected for double 
lenses are an order of magnitude smaller than in quad lenses, so it comes as no surprise 
that there is no convincing determination of a relative proper motion in the 0957+561 lens 
(Campbell 1995, Campbell et al. 1994). 

The proper motions are dominated by the proper motions of the lens galaxies, because 
most estimates of the rms peculiar velocities of galaxies are significantly larger than our 
motion relative to the CMB frame. If the rms peculiar velocity today is Vrms,o and we 
assume it evolves with redshift as v,rms{z) = ^^rms,o/(l + ^Y^'^^ then the time scale for a 2-(T 
detection in a typical quad lens is 

Dos 



T = 10 



600 km 



(^0 



4:g{zs) Gpc 



yrs, (31) 



10/ias 

assuming that the rms uncertainty in the measurement of one coordinate of an image 
position is ag = 10 //as (optimistic), but that only two measurements are made, and using 
only one image pair instead of all six independent measurements (pessimistic). Moreover, 
our statistical model uses an underestimate of the rms magnification of the proper motions. 

Given an ensemble of quad lenses with proper motions we can estimate the proper 
motion dipole due to our motion relative to the extragalactic inertial frame. In any one 
system, the typical magnitude of the proper motion dipole is smaller than the peculiar 
dipole, but it is detectable when averaged over an ensemble of lenses. Given ten lenses and 
an rms uncertainty of ag = 10 //as, we estimate that the direction of the dipole can be 
determined to 20° and the Hubble constant can be determined to AHq/Hq = 0.34 after ten 
years (1-cr limits). After ten years the uncertainties are dominated by the motions of the 
lens and the source rather than measurement error, so precision measurements depend on 
finding larger numbers of lenses. 
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Barrow, Juszkiewicz, and Sonoda (1985) put limits on the dimensionless quantity 
lo/Hq by using constraints from the quadrupole moment in the CMB multipole expansion. 
These hmits are still relevant in view of the current COBE results (Bennett et al. 1994) 
of lo/Hq < 2 X 10~^ (10""*) for a flat (open) universe. In turn, these numbers permit 
a maximum rotation rate of 10~^'*rad yr~^ ~ 0.02 //as yr~^. The rotation limits of the 
"universe" (i.e. our horizon) from CMB quadrupole [u ^ 0.02/ias yr~^) are stronger than 
the possible limits from the current proposed analysis. Rotation, however, can also be 
contributed from a broken anisotropy on smaller scales than the horizon. In the latter case 
the current method can put limits of the order to < 10~^^ radyr"^. 

The systematic uncertainties in this approach are dominated by the poor limits on 
the internal motions of VLBI sources that are not obvious superluminals. If these motions 
exceed ~ 2000 km s~^, they will dominate the proper motions on short time scales. Internal 
motions and microlensing do not affect the longer term proper motion of the lens, because 
both processes only introduce extra noise in addition to the monotonic proper motions 
produced by the velocity of the observer, lens, and source galaxies. Existing estimates for 
proper motions of extragalactic VLBI sources are consistent with zero given their systematic 
errors (e.g. Marcaide et al. (1994), Campbell (1995), Campbell et al. (1994)). If, however, 
typical sources had apparent motions even 10% of the values seen in 1038+528 A/B (31 ±22 
//as yr~^) or 0957+561 (~ 25 //as yr~^) under the assumption that the motion is real, then 
the VLBI quad lenses should show astounding proper motions (~ 300 //as yr~^). While we 
view this as unlikely, such proper motions would be extraordinarily useful for constraining 
the lens geometry and probing AGN physics. 

We would like to thank L. Greenhill, J.N. Hewitt, J.M. Moran, and I.I. Shapiro for 
discussions, and J. Miralda-Escude and U. Seljak for pointing out an error. This work was 
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